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Abstract 

We present the threshold resummation for the cross section for electroproduction of 
heavy quarks. We work to next-to- leading logarithmic accuracy, and in single-particle 
inclusive kinematics. We provide next-to-leading and next-to-next-to-leading order ex- 
pansions of our resummed formula, and examine numerically the quality of these finite 
order approximations. For the case of charm we study their impact on the structure func- 
tion Fi and its differential distribution with respect to the charm transverse momentum. 



1 Introduction 



The production of heavy quarks in deep-inelastic scattering is a reaction of great interest because 
it sheds light on a number of fundamental issues in the QCD-improved parton model. First, 
at moderate momentum transfer Q, it is a direct probe of the gluon density in the proton. 
Second, it allows the exploration of the relevant degrees of freedom, as a function of Q, for 
describing this process, the options comprising the treatment of the heavy quark as a quantum 
field external to the proton, via an "effective" parton density, or as a combination thereof. 
Finally, detailed studies of perturbative QCD dynamics are possible using the presence of at 
least two hard scales: Q and the heavy quark mass m. 

Charm quarks have been identified in deep-inelastic electron proton collisions by the EMC 
collaboration |l| at small Q and low final state invariant mass (large x), and more recently and 
in much greater numbers at high Q and medium to small x by the ZEUS and HI collaborations 
at HERA ||. Considerably more charm (and bottom) data are anticipated at HERA. 

At the theoretical level the reaction has already been studied extensively. In the framework 
where the heavy quark is not treated as a parton, leading order (LO) || and next-to-leading 
order (NLO) [| calculat ions of the inclusive structure functions exist. Moreover, LO f| and 
NLO calculations of fully differential distributions have been performed in recent years. So 
far theory agrees reasonably well with the HERA data. We shall not concern ourselves here 
with descriptions in which the heavy quark is (partly) treated as a parton j7|. 

At the HERA center of mass energy of 314 GeV, one naively expects the production of 
a heavy quark pair with a minimum invariant mass of a couple of GeV to be insensitive to 
threshold effects, however a closer examination || reveals that this is not true: the density 
of gluons, which controls the heavy quark electroproduction rate, is large at small momentum 
fractions of the proton, favoring partonic processes where the initial gluon has only somewhat 
more energy than needed to produce the final state of interest. 

As we shall show, for the inclusive cross section and for x values larger than about 0.01, the 
higher order QCD corrections for this reaction are in fact dominated by singular distributions 
(a"/n!) \hi 2n ~ (w)/w] + at order n. The dimensionless weight to is a function of the momenta of 
the external partons, chosen in accordance with the kinematics @, and vanishes at threshold. 
Our main purpose in this paper is to resum these singular functions to all orders in perturbation 
theory to next-to-leading logarithmic (NLL) accuracy, and, moreover, in single-particle inclusive 
(1PI) kinematics. The required technology has recently been developed in Refs. [[]|0],|n].|I^. 

We also provide analytic results for the resummation and finite order expansion for pair 
invariant mass kinematics (appendix B). All our analytic results are valid for either charm or 
bottom production, but we concentrate our numerical studies on the case of single inclusive 
charm production. 

The interest in this analysis is first of all intrinsic, and lies in studying the quality of the 
next-to-leading logarithmic resummation for single-particle inclusive observables. Second, while 
the exactly computed 0(a s ) corrections to the inclusive structure function ,F| harm [Q are, for 
experimentally accessible x values, on average about 30% - 40% [fL4|], i.e. non-negligible but 
not too alarming, the size of even higher order corrections bears examination. Our results, in 
particular our estimates for the unknown exact next-to-next-to-leading order (NNLO) O(a^) 
corrections, should help establish the theoretical error on observables for the reaction under 
study, and on the gluon density extracted from charm quark production at HERA. Finally, our 
results might contribute to a future global analysis for resummed parton densities. 

We have organized our paper as follows. Section 2 contains the derivation of the resummed 
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formula for the single-particle inclusive differential cross section for heavy quark electroproduc- 
tion. In section 3 we study its finite order expansions both analytically and numerically. We 
examine both partonic quantities (coefficient functions), as well as hadronic observables, viz. 
the inclusive structure function i^ 2 charm and its differential distribution in charm transverse mo- 
mentum. Section 4 contains our conclusions. In appendix A we present some useful formulae 
related to Laplace transformations, and in appendix B we perform the threshold resumma- 
tion and NNLO expansion for the heavy quark electroproduction cross section in pair-inclusive 
kinematics. 



2 Resummed differential cross section 

We begin with the definition of the exact single-particle inclusive kinematics. We study the 
electron (e) - proton (P) reaction 

e(l) + P(p) - e(l -q) + Q( Pl ) + X[Q](p' 2 ) , (1) 

where X[Q] denotes any allowed hadronic final state containing at least the heavy antiquark, 
Q(pi) is a heavy quark, and g M the momentum transfer of the leptonic to the hadronic sector of 
the scattering. After integrating over the azimuthal angle between the lepton scattering plane 
and the tagged heavy quark plane, neglecting Z-boson exchange (we assume \q 2 \ <C Mf ) and 
summing over X, the cross section of the reaction ([[]) may be written as [15 



dxdQ 2 dT 1 dU 1 xQ 4 



(l + (l-y) 2 ) 



rf 2 F 2>P (x, Q 2 , m 2 , 71, Ui) 2 d 2 F L , P (x, Q 2 , m 2 , 71, U x ) 



y 



dTi dUi dT\ dU\ 

(2) 

The functions d 2 F k pjdT\dU\, k — 2,L are the double-differential deep-inelastic heavy quark 
structure functions and a is the fine structure constant. The kinematic variables Q 2 ,x,y are 
defined by 

Q 2 = -q 2 >0 , x = ^- , y = P -^. (3) 

2p ■ q p ■ I 

We also define the overall invariants 

S = (p + q)* = S'-Q 2 , T 1 = (p- J 9 1 ) 2 -m 2 , Ul = (q- Pl ) 2 -m 2 , (4) 

and 

S 4 = S + 71 + U x + Q 2 ■ (5) 

The structure functions d 2 F k>P /dTi dU\ describe the strong interaction part of the reaction ([[J). 
They enjoy the factorization 

d 2 Fk,p(x, S4, Ti, Ui, Q 2 , m 2 ) 1 ^ r dz . 2 rx s± t\ U\ Q 2 m 2 



f dz 2 ( x s 4 ti Mi Li" mf . , 

2^ / —9i/p{z,H ) uj k A-, — , — , — , — , — ,a s {ii) 
.= J z K z u z u 2 u z u z u l 



dTidUi S' 2 iJ^qj J_ z ' ' ^z' /J 2 ' /i 2 ' /i 2 ' [i 2 ' [l 2 

(6) 

where the sum is over all massless parton flavors, and (pi/p(z, n 2 ) is the parton distribution 
function (PDF) for flavor i in the proton and z the momentum fraction of parton % in the 
proton with z~ = —U\/(S' + 71). The dimensionless functions u>k,i describe the underlying 
hard parton scattering processes and depend on the partonic invariants s',ti,ui and S4, to be 
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defined in Eq. (|8|). The factorization scale is denoted by /i and in this paper is always taken 
equal to the renormalization scale. 

As we are operating near threshold we will assume a fixed number of light flavors in the 
evolution of the 4>i/p and the strong coupling a s . A further simplification we adopt is neglecting 
the contributions from quarks and antiquarks in the sum over flavors in Eq. This is well- 
justified as their contribution at NLO was found to be a mere 5% {4,14]. We therefore only 
consider the gluon-initiated partonic subprocess 



7*0?) + g(k) - 

where k = zp. The partonic invariants are 



Q( P i) + X'[Q}(p' 2 ) 



s = (k + q) 2 =s' + q 2 



ti = (k — pi) 2 — m 2 



Ui 



(q — Pi) 2 ~ m<2 



(7) 



(8) 



Note that t\ = zT±, s' = z(S + Q 2 ) and U\ = U\. The invariant s 4 = M\, — m? measures the 
inelasticity of the reaction (0) and is given by 



s' + ti + U\ 



s 4 • 



(9) 



Our goal is to resum those higher order contributions to u)k, g that contain 1PI plus-distributions. 
The latter are defined by 



"ln'(s 4 /m 2 )~ 


= lim < 


s 4 





ln / (s 4 /m 2 ) 



«4 



(10) 



These distributions are the result of imperfect cancellations between soft and virtual contri- 



butions to the cross section. At order 0(a 



i+2\ 



0, 1, . . . the leading logarithmic (LL) 



corrections correspond to I — 2i + 1, the NLL ones to I = 2i etc. 

In pursuing our goal, we follow the general principles of Ref. 0, and the methods of 
Refs. |l^JTT| , [T2|JT6| . We restrict ourselves from here onwards to the structure function F^p. 
The resummation of F^ p demands special attention [Hj. Moreover, F 2) p constitutes the bulk 
of the cross section in Eq. (|). 

By replacing the incoming proton in Eq. (0) with an incoming gluon, one may compute the 
hard part cl> 2j9 in infrared-regulated perturbation theory. The resummation of u>2, g rests upon 
the simultaneous factorization of the dynamics and the kinematics of the observable in the 
threshold region of phase space. This factorization is pictured in Fig. [I]. The figure corresponds 
to the partonic process (0) and represents a general partonic configuration that produces large 
corrections, i.e. corrections containing the singular functions of Eq. (|T(i|). 

The figure shows the factorization of the partonic cross section into various functions, each 
organizing the large corrections corresponding to a particular region of phase space. Such 
factorizations have been discussed earlier for deep-inelastic scattering, Drell-Yan [[L8[], heavy 
quark |T6|JT9|] , dijet [T(J and for general single-particle/jet inclusive cross sections fT^fl . The 
function ip g / g contains the full dynamics of partons moving collinearly to the incoming gluon g. 
It includes all leading and some next-to-leading enhancements. The function S(ks) summarizes 
the dynamics of soft gluons that are not collinear to g, and includes all remaining next-to- 
leading contributions. The function B.2,g = h 2 g h<2,g incorporates the effects of off-shell partons, 
and contains no singular functions. There are large next-to-leading corrections associated with 
the outgoing heavy quarks. Were we to treat the heavy quarks as massless, each heavy quark 
would be assigned its own jet function, see e.g. Refs. ||i0yi3||l"3| . In our case the heavy quark 
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mass prevents collinear singularities, so that all singular functions arising from the final state 
heavy quarks are due to soft gluons. Hence, near threshold, all the singular behavior associated 



with the heavy quarks may be included in the soft function [I6|. We note that S(ks) is simply 
a function, and not a matrix in a space of color tensors, in contrast to heavy quark |16| or jet 
production p~0|JTT[| in hadronic collisions. 



The kinematics of reaction (|7|) decomposes in a manner that corresponds exactly to the 
factorization of Fig. [TJ. Momentum conservation at the parton level means 

qJ rZ p = p 1 J r p 2 J r k s , (11) 

where, in view of the above discussion, pi and p 2 may be interpreted as the on-shell momenta 
of the heavy quark and anti-quark respectively. Squaring and dividing by m 2 , we have 

Si M _ x( 2p- q' -2p-p^ 2p 2 -k s 



m 2 ^ m 2 ' m 2 



(l-z)(^) + ^= Wl (^) +Ws , (12) 
v in 2 ' m 2 v m 2 ' 

?2 



where we have dropped terms of order S\. Notice that the kinematics is in fact specified by 
the dimensionless vector defined as Q 1 = p^/ni- In Eq. (|T^) we have identified the overall 
single-particle inclusive weight w, mentioned in the introduction, with S^/m 2 . At fixed 64 the 




Figure 1: Factorization of heavy quark electroproduction near threshold. The double lines denote 
eikonal propagators. 



infrared-regulated, differential partonic structure function factorizes as 
, 2 d 2 F 2 , g (x, S^tuuu Q 2 ,m 2 ) 2 2 r 

s df~du = H2 >o( tl > Ul >Q > m ) J dwidws (13) 

x£ (-^ ~ w i(-^p) ~ w s) ^9/ a iw u pX,ri) S(w s ,f3 h (,n) , 

where the kinematics relation (|l^) is implemented in the ^-function. The $ are the four- 
velocities of the particles in the Born approximation to Eq. (0), i.e. with X[Q](p 2 ) replaced by 
Q(p 2 )- The various functions are computed in n ■ A = gauge. 

Further manipulations are carried out most conveniently in terms of Laplace moments, 
defined by 

00 

f(N) = fdwe- Nw f(w). (14) 
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The upper limit of this integral is not so important, and may also be put at 1, where \nw = 0. 

As both the gauge and kinematics vector are timelike, we may choose [] = in which 
case the 1PI ipi/i are equal to the center-of-mass densities (for which ^ = n M = 5^°) [O]. 

Replacing the incoming proton with an incoming gluon in Eq. 
Eq. (|13|), we derive 



and comparing with 




«i Q 2 m 2 




m 



x 



jLf /i z [I* fJL* J 



C/aO) 



4>g/g{^ '(-Ux/m 2 ) ,{Jt) 



s{ 111 



(15) 



ensure that the N- 



Given the factorization in Eqs. ( pH) and (0), the arguments of 
dependence in each of the functions of Eq. (fL5|) exponentiates. 

We next discuss each function in Eq. ([IB]), or (|T|), in turn, starting with the ip g / g wavefunc- 
tion. In analogy with the center-of-mass definitions [18,10 it may be defined as an operator 
matrix element f\ 

on 

c 



ll>g/g(w,P,(,n) 



1 



P 



2{N 2 - l)^(vp) 2 
x (g( P )\F^(yO 



dye- iy{1 - w)p< 
Lv (0)\g(p)) n .A= , 



(16) 



where w = 1 — z and is a fixed lightlike four- vector in the opposite direction as p^, such 
that v ■ p is of order one. The first factor refers to a spin and color average. The states are 
normalized such that (0|v4 J -(0)|g(£>)) = e^Qo). Expression (|T6D as a whole is normalized such 
that 



S(w) 



(17) 



We have calculated the order a s corrections to this operator matrix element in general axial 
gauge n • A = 0. The 0(a s ) corrections to the operator matrix element in ( |T6"D are in d = 4 — 2e 
dimensions, and to next-to-leading logarithmic accuracy 



/(!)/- p 



c 




2 ln(tu) 



w 



- 1 + In(2i/ a ) + 0(e) 



(18) 



where z/ g = ((3 g ■ n) /\n\ , with n M chosen equal to and the four-velocity of the incoming 

gluon, defined by p>* = f3^Js/2. Note that \n(2u g ) + ln(s/// 2 ) = ln(4(p • C) 2 /V)- We have 
abbreviated 1/e = 1/e — 7s + ln47r. We have verified that the steps for the resummation of the 
Sudakov double logarithms ln 2 (iV) in this function trace exactly those for the center-of-mass 
density introduced for the Drell-Yan process in Ref. [JTSj . 

We choose for <p g / g the MS density, which does not have Sudakov double logarithms. To the 
same accuracy as ip g / g in Eq. fllSD , it is given by 



<frg/g (W,« s (/i), 



6{w) + 



Oi s (jj) 



71 



(t) 


V 




.w. 



(19) 



■"■Note that for lightlikc kinematics, such as in direct photon production, it may be convenient to keep both 
vectors separate frL2| . 

2 For the equivalent definition for (anti)quarks, see Ref. J12[. 
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The resummation of the singular contributions to this function is done via the Altarelli-Parisi 
equation [pOfl , see also BUI- We only need this result in combination with the resummed ip g / g 
as the ratio ipg/ g /4> g /g in Eq. (|15|). This ratio is finite, as can readily be verified at one-loop 
from Eqs. (|18p and flf|). 

The scale dependence of these functions may be found as follows [[UJ • The scale dependence 
of ipg/g is governed by 



d ( p ■ ( \ 



(20) 



The anomalous dimension on the right does not depend on the moment N, as the only ultraviolet 
divergences are due to gluon wavefunction renormalization. Hence, 7^ = 2 / y g , where 7 9 is the 
anomalous dimension of the gluon field in axial gauge. 

The anomalous dimension of the MS density does depend on N, via the Altarelli-Parisi 
equation 

V-j^^(f) g/ g(N,a s (iJ,),e) = 2-f g / g (N,a 8 (fi)) . (21) 

We next discuss the soft function S(N,Q. It summarizes the singular contributions arising 
from soft gluons that are not collinear to the incoming jet, and is infrared finite. It depends 
on the kinematics of the hard scattering and contributes only at next-to-leading logarithm. 
Threshold resummation for soft functions is treated in detail in Refs. [0,01; a brief sketch will 



suffice here. The factorization in Eq. (|13D introduces ultraviolet divergences, distributed in such 
a way between the hard function H and the soft function S, that they cancel in the product 
P,|T0f. The function S may be defined as the matrix element of a composite operator that 
connects Wilson lines in the directions of the external partons. Its extra ultraviolet divergences 
are cancelled by the renormalization of this operator. As a result, S obeys the renormalization 
group equation 



' d 
li— + P(a s 



d 

■ ■ 

da. 



S = -(2Re{T s {a s {fi))})S, 



(22) 



then resums the ln(iV) terms in S. In the reaction under study the anomalous dimension is a 
lxl matrix in color space, so we may solve Eq. (22) directly: 



exp 



/ 42Re{r 5 (a s (//))} 



(23) 



The resummed hard part may now be written in moment space as 
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The first exponent summarizes the iV-dependence of the wave function ratio ipg/g/fig/g an d 
is, in moment space, the same as for heavy quark and dijet production pj| , |lT 



E (g) ( N u, 



m 



dw 



-N u w\ 



w 



+ 2^(9) M wm )} 



(25) 



"GO = 2^-|l-ln(^— _ 
where N u = N(—ui/m 2 ), and is given by the expression 



^)M = c a [- + -k[- 



(26) 



(27) 



with K = Ca (67/18 — 7r 2 /6) — 5/972 j [Q and the number of quark flavors. 

The second exponent controls the factorization scale dependence of the ratio ipg/g/<t> g / g 
through the anomalous dimensions 7 S and J g / g - In axial gauge we find them to be 



"ig/g (N u ,a s (n)) 



a s (n) 



7T 



a s {jj) 



71 



(C A ln(N u ) - 6 2 ) + 0(l/N) . 



(28) 
(29) 



where b 2 = (11C A - 2n f )/12. 

Finally, the 0(a s ) expression for the soft anomalous dimension T$ can be inferred from the 
UV divergences of the eikonal Feynman graphs in Fig. 0. We obtain the result 




Figure 2: One-loop corrections to the soft function S for heavy-quark production in photon-gluon 
fusion. The double lines denote eikonal propagators. The lines labelled k corresponds to a gluonic 
eikonal line, and those labelled p± and p 2 denote quark and antiquark eikonal lines, respectively. 



7T 



C f 



1 - In 



rrr 




(30) 



where L@ is given by 



1-2^7- f^WV. 4 , ^=713^, (31) 



(3 \ + 
and, for single-particle inclusive kinematics, 

m m 2 



(32) 



At the conclusion of this section a few remarks are in order. First, the integral in the 
exponent in Eq. (|25|) can only be interpreted in a formal sense, as some prescription should 
be implemented to avoid integration over the Landau pole in the running coupling. We do 
not address such renormalon ambiguities here. In this paper we only employ the resummed 
expressions as generating functionals of approximate perturbation theory. 



Second, although our main focus is on 1PI kinematics, it is straightforward [|T^] to derive 
from Eq. ([24]) the resummed expression for pair-invariant mass kinematics, i.e. the resummed 
cross section for the process in Eqs. (ffl) and (|7|), differential in the QQ invariant mass. We 
present this derivation in appendix B. 



3 Finite order results 

In this section we expand our resummed cross section to one and two loop order so as to 
compare with exact NLO calculations and provide NNLO approximations. 

We begin by deriving NLL analytic formulae for the single-particle inclusive, partonic hard 
part 0J2, g - Subsequently we study their behavior in the hadronic inclusive structure function 
F 2 charm , and in the differential distribution dF 2 charm /<ip T . 



3.1 Partonic results at NLO and NNLO 

The Born level hard part (see e.g. for the process 7* + g — > Q + Q is 

uji%',t hUl ) = 5(s' + t l +u l )a^(s'M,u l ) ) 

with cr 2 3 ° rn given by 



(33) 



a 



Born / „/ 



2.g 



Os_ 2 N C C F 2 

27i €<1 iV c 2 - 1 W 



ti u\ . m 2 s' 

- + -+ 4 T 



/ 9 A 9 9 

s q „ q n m q 



,'2 



m 2 s' s 
t 1 u 1/ 

J2 



2 J 



till! 



2 - 



+ 



S' \ t\ U\ 



s 
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(34) 



We now derive the NLO soft gluons corrections to u>2, g (s', ti,ui) by expanding the resummed 
hard part in Eq. ([24]) to one loop, using the explicit expressions for the various functions in 
Eq. (|D, given in Eqs. ©, (|l|) and We find 



a 2 B ° r V,t 1)Ml ,, 



(35) 
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where 



+ 



W\=Si/{ — Ml),£^0 



w s 



a s (fi) 



2C f 



ln(s 4 /m 2 

Si 



+ 



1 

Ls 4 



2Re{r s } 



C A {ln(^)+ReL B -\n{ iL 



(36) 



■ws=S4,/m 2 



-2 C F (ReL^ + lH + <5(s 4 ) CU In f In f ^4 



(37) 



with fi the MS-mass factorization scale and Lp given in Eq. (pip. The plus-distributions in s 4 
have been defined in Eq. (p~0|) . To next-to- leading logarithm, the result in Eq. fl3~7|) agrees with 
the exact a s corrections of Ref. ||]. This holds for the contributions both independent and 
dependent on the factorization scale, i.e. terms including powers of ln(/x/m). 

The scale-dependent logarithms ln(/x/m) are explicitly generated by the second exponent in 
Eq. (p^D, involving the anomalous dimensions 7 g and 7 9 / 9 given in Eqs. ( pB| ) and (P5|). In deriving 
7 9 / g in Eq. (|29"D we have carefully determined all terms which are constant in moment space, 
including the ln(— u\/m 2 ) logarithm, absorbed in N u . In this way we are able to determine all 
factorization scale dependent terms to NLL accuracy. 

The NNLO corrections may be found in a manner exactly analogous to what was done for 
Drell-Yan in Ref. |23 and for Higgs production in Ref. |]2T|. We expand Eq. (|24] ) to second 
order in a s in moment space. With the table of Laplace transforms in appendix A for the 
singular distributions in Eq. (|10"D, we then find the 0{a 2 s ) threshold-enhanced corrections in 
momentum space. Our result is 



(38) 



with 



7T 



2C\ 



In 3 (s 4 /m 2 ) 



(39) 



+ 



r ln 2 (s 4 /m 2 ) 



s 4 



In 



"i 



+ ReL 



hi 



-2C A (6 2 + 3C F (ReL /3 + l)) 



+ 



ln(s 4 /m 2 

Si 



b fe <- 21n (£)- 2R ^ + 21n (9 



+2C A (b 2 + 2C F (ReL (3 + l)) 



+ 



1 

S 4 



ln 2 (^H-C 2 ln 



-Ml 



and b 2 = (11C A - 2n/)/12. 
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3.2 Gluon coefficient functions 

The inclusive structure function F^p rra is obtained from Eq. (|6]) by integrating over T\ and U\. 
It can be expanded in coefficient functions c 2g as follows (dropping the P subscript) 

2/02 \ oo k ,,2 



Ff-(*,Q 2 ,m 2 ) = a fl el f dzJ> g/P (z^ 2 ) t(^<»)) k Y,$ 9 l \v,tW— 2 , (40) 

k=0 1=0 



where a = (Q 2 + 4m 2 ) /Q 2 and we recall that contributions from light initial state quarks are 
neglected. The function <j) g /p(z, fi 2 ) denotes the gluon PDF and the functions c 2 J depend on 
the scaling variables 

s i, e = ^. (4i) 



4m 2 m 2 

The variable rj is a direct measure of the distance to the partonic threshold. The inclusive 
coefficient functions are obtained from Eqs. ( |5?| ) and ( |39"D by 



s'(l-/3)/2 



where the double differential coefficient functions are in turn related to the hard part lu 2j9 of 
Eq. (|) by 



^ /2 ^4yV,*i,'"i) , i V 2 _ 4tt 2 m 2 (fc) 

^ S dti<fc 4 ™ 2 " a s e 2 (4vra^ Q 2 <J ^ 2 ^ ^ 

Also, in Eq. © we abbreviated [|] 



(43) 

«l=s'+tl— S4 



We begin our numerical studies with the gluon coefficient functions c 2g (rj, £), i.e. those 
that are not accompanied by scale-dependent logarithms. We recall that all our results are 
derived in the MS scheme. The functions c^f^ and c 2 ]g ^ are known exactly, see Ref. fij. The 
one-loop expansion of our resummed hard part provides an approximation to the exactly known 
c 2>g . To judge the benefits of resumming to next- to -leading logarithmic as compared to leading 
logarithmic accuracy, we can distill also LL expressions, by keeping only the [ln(s4/m 2 )/s 4 ] + 
and [ln 3 (s 4 /m 2 )/s 4 ] + terms in Eqs. (|37| ) and ([39]) respectively. 

In Figs. |a and [|a we plot the functions c 2g { (77, £), k = 0, 1 versus 77, for two values of 
£ = Q 2 /m 2 . For the coefficient functions only the ratio £ matters, however, we chose those 
values to correspond to Q 2 = 0.1 and 10 GeV 2 for a charm mass [] of m = 1.5 GeV. These 
figures reveal that, at one loop, the LL accuracy provides a good approximation for very small 
77, however, with significant deviations from the exact result for larger 77, i.e. as the distance 

3 In the next subsection, where we study the hadronic structure function, we use the value m = 1.6 GeV. 
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Figure 3: (a): The ^-dependence of the coefficient functions (r],£), k = 0, 1 for Q 2 = 0.1 GeV 2 

with m = 1.5 GeV. Plotted are the exact results for c^ g °\ k = 0, 1 (solid lines), the LL approximation 

to C2g^ (dotted line) and the NLL approximation to c^'^ (dashed line), (b): The ^-dependence of the 

coefficient function Cgp (rj, £) for Q 2 = 0.1 GeV 2 with m = 1.5 GeV. Plotted are the LL approximation 
(dotted line) and the NLL approximation (dashed line). 
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Figure 4: (a): The ^-dependence of the coefficient functions cj£' ° (?7,£), k = 0,1 for Q 2 = 10 GeV 2 
with m = 1.5 GeV. The notation is the same as in Fig. |3|a. (b): The ry-dependence of the coefficient 
function c^L (rf, £) for Q 2 = 10 GeV 2 with m = 1.5 GeV. The notation is the same as in Fig. |3|b. 
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Figure 5: (a): The NLL approximation to the coefficient function c 2g for Q 2 = 10GeV 2 and 

m = 1.5 GeV with restrictions to the small ry-region. Plotted are the unmodified NLL result (solid 
line), the NLL result with a factor 9{m 2 — S4) included (dotted line) and the NLL result multiplied 
with a damping factor + r] (dashed line), (b): The ^-dependence of the NLL approximation to 

the coefficient function c^'g (v>Q fo r Q 2 = 10 GeV 2 and m = 1.5 GeV. The notation is the same as 
in Fig. ||a. 



from threshold increases. On the other hand, the NLL approximation is excellent over a much 
wider range in 77, up to values of about 10. 

We also show c^ g °\ obtained from Eqs. (|39"D and fl43|) , in LL and NLL approximation in 
Figs. [3]b and [|b for the same values of m and Q 2 as before. We observe more structure than 
in the c^ g °^ curves. Figs. [|b and (|b show some sizable deviations from zero at large 77, which 
leads us to consider the following issue. Kinematically, larger 77 values allow w = s 4 /m 2 3> 1, 
whereas, in contrast, the threshold region is denned by w ~ 0. Therefore we should consider 
the inclusion of a factor 6(m 2 — S4) in the integral ( fi2| ) in order to suppress spurious large 
corrections in the large 77 region. 

In Figs. |5]a and |5|b we show the effect of the veto 6(m 2 — s 4 ) on the NLL coefficient functions 
for k = 1, 2 for values of Q 2 = 10 GeV 2 and m = 1.5 GeV. Also shown is the effect of 
a damping factor + rj p^| , which can be included outside the integral in Eq. (f4^ ) and 

is therefore easier to implement. We see that the effect of the veto is well- mimicked by the 
damping factor, therefore we shall use the latter, rather than the veto, where needed in what 
follows. Of course, for smaller values of 77 the damping factor leads to a slight underestimation 
of the coefficient functions compared to the veto, but these effects are negligible, as Fig. |] 
shows. 



Let us now investigate those coefficient functions c^g (77, £), / > 1 in Eq. (^OJ) that determine 
the dependence of _p^ harm on the mass factorization scale \i. As mentioned earlier, this is of 
some relevance since the exact NLO _F£ harm exhibits considerably more /1 sensitivity at large x 
(near threshold), than at small x, cf. Ref. ||. The natural question arises whether and, if so 
how much, our approximate NNLO results might improve the situation. 

Explicitly, the coefficient functions under consideration up to NNLO are c^P which is 

known exactly 0] and the previously unknown functions c^'p and . We can infer the exact 
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Figure 6: (a): The 77-dependence of the coefficient function c 2 ^ g (rj,0 f° r Q 2 = 0.1 GeV 2 with 

m = 1.5 GeV. Plotted are the exact result (solid line), the LL approximation (dotted line) and the NLL 

(2 1) 

approximation (dashed line), (b): The //-dependence of the coefficient functions c 2g (rj, £), I = 1, 2 for 
Q 2 = 0.1 GeV 2 with m = 1.5 GeV. Plotted are the exact results (solid lines), the LL approximations 
(dotted lines) and the NLL approximations (dashed lines). 
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Figure 7: (a): The ^-dependence of the coefficient function c 2 g (77, £) for Q 2 = 10 GeV 2 with 
m = 1.5 GeV. The notation is the same as in Fig. ||a. (b): The ^-dependence of the coefficient 
functions c£'g\rj,€), 1 = 1,2 for Q 2 = 0.1 GeV 2 with m = 1.5 GeV. The notation is the same as m 
Fig. |b. 
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Figure 8: (a): The NLL approximation to the coefficient function c 2g (t),£) for Q 2 = 10GeV 2 and 
m = 1.5 GeV with restrictions to the small ^-region. The notation is the same as in Fig. [|a (b): The 
NLL approximation to the coefficient function c^ g l \r},£), ^ = 1,2 for Q 2 = 10 GeV 2 and m = 1.5 GeV 
with restrictions to the small r/-region. The notation is the same as in Fig. [|a. 



results for these functions from renormalization group arguments. We find 
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(45) 
(46) 
(47) 



with /^-function coefficients b 2 = (11C A - 2n f )/12, b 3 = (34(7 2 - QC F n f - 10C A n f )/48 and the 
one- and two- loop gluon-gluon splitting functions Pgg\ Pgg i c f- R^f- BUI- The convolutions 
involving a coefficient function c 2 ^ are defined as 



(48) 



where a = (Q 2 + 4m 2 )/Q 2 and 77(2;) is derived from Eq. (fllD, ^(x) = £/4(l/x — 1) — 1. Finally, 
the function in Eq. (g7|) is given by the standard convolution of two splitting functions 



P(0) 
99 1 



1 

dxx | dx 2 5(x - x 1 x 2 )P^(x l )P^(x 2 ) . 



(49) 



With Eqs. (]4"5|)-(f4"T|) in hand, we are able to check the quality of our NLL approximation 
at NLO and NNLO against exact answers. Again, to fully appreciate the improvement of NLL 
resummation, we truncate our approximation to LL accuracy, keeping only the ln(/x/m)[l/s4] + 

4 Eq. ( F^[ ) agrees with the result for c^'^ in Ref. [Q. 



15 




10 " 3 10 " 2 10 _1 1 io io 2 io 3 10 " 3 10 " 2 io" 1 1 io io 2 io 3 



r) = s/(4m 2 ) - 1 r| = s/(4m 2 ) - 1 

Figure 9: (a): The CTEQ4M gluon density 4> g /p{x/z, /x 2 ) as a function of rj for m = 1.6 GeV, 
Q 2 = 10 GeV and x = 0.1. Plotted are the scale choice /i = m (solid line), \i = 2m (dashed line) and 
fj, = \/Q 2 + 4m 2 (dotted line), (b): Same as Fig. ||a for x = 0.01. 



term in Eq. fl37|) and the terms In(/i/m)[ln 2 (s4/m 2 )/s4] + and ln 2 (/i/m) [In(s4/m 2 )/s4] + in 
Eq. ([59|), respectively. 

In Figs. |B|a and 0a we plot the function c^' g (tj,£) versus 77, for two values of £ = Q 2 /m 2 , 
while Figs. |6|b and [F]b display the functions c^' g l \ I — 1,2. The evaluation of Eq. (|46|) has 



been performed with the help of the parametrization for c^'^ from Ref. [p6 |. From Figs 



and 0, we see again that the NLL approximations are superior to the LL ones. In Figs. [5]a 
and [8|b we show the effect of the veto 9{m 2 — S4) in the integral (|42|) and the damping factor 
+ 77 on c^ 1 ^//, £) and c^, (77, Z = 1, 2 for values of £ corresponding to Q 2 = 10 GeV 2 
and m = 1.5 GeV. In the following sections we shall for the sake of uniformity include a 
damping factor for the ^^(r],^), I = 1, 2 and (^^(t], £), even though the undamped NLL c^ 1,1 ' 



.'7 



approximates the exact somewhat better, as a comparison of Figs. 0a and [|a shows. We 
verified that the difference is negligible at the hadronic level. 

3.3 Inclusive structure function F 2 charm 

Having obtained encouraging results for the quality of the partonic approximations Eqs. ( |3T| ) 
and (|39D at the inclusive level, we shall now assess their effect at the hadronic level. 

Throughout we use the CTEQ4M f27j gluon density. For NLO plots we use the two-loop 



expression for a s with n/ = 4, A = 0.298 GeV. Respectively, for NNLO plots, the three-loop 



expression |28[ for a s and an adjusted A = 0.265 GeV. 

In Fig. |9]a and []b we plot the gluon density <f> g /p(x/z, /x 2 ) as a function of rj on the same 

scale as the coefficient functions c^'g (t},£) for x = 0.1 and 0.01 and three different choices 
of the factorization scale (/i = m, 2m, \JQ 2 + Am 2 ). From now on we use a charm mass of 
m = 1.6 GeV and take Q 2 = 10 GeV. Our heavy quark mass is always a pole mass, as in 
Ref. ||]. The figure shows that the gluon density indeed provides support for the coefficient 



functions in Eq. (EO) in the threshold region for x > 0.01, although the figure also shows the 



extent of the support to be somewhat scale sensitive. 

A more revealing way to examine the support of the gluon density in Eq. fl40D is to plot the 
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Figure 10: (a): F 2 charm (x, Q 2 , z max ) as a function of z max at NLO with the CTEQ4M gluon PDF, 
x = 0.1, m = 1.6 GeV, Q 2 = lOGeV and /x = m (upper three curves), /U = \/Q 2 + 4m 2 (lower three 
curves). Plotted are: The exact results (solid lines), the LL approximations (dotted lines) and the 
NLL approximations with the damping factor 1/^/1 + r\ on the scale dependent terms only (dashed 
lines), (b): Same as Fig. |lO| a for x = 0.01. The lower three curves, \i = \JQ 2 + 4m 2 , have been scaled 
down by a factor of 2. 



integral as a function of z max , 

( \ p 2f)2 2n }? x oo k 2 

F 2 ch -(x,Q 2 ,m 2 ,z max ) ~ a f / ^0 s/ p(z,/i 2 ) E( 47r «.(A*))*E^(^ 0^-^.(50) 

ax re— u t— u 

Varying the theoretical cut-off z max allows us see where the integral (pf^) acquires most of its 
value. The physical structure function corresponds to z max = 1. 

In Figs. [10| a and |10| b we show the results for the numerical evaluation of Eq. ( |50"D and 
compare to the exact NLO calculation of Ref. [|J for two values of x — 0.1,0.01 and two 
choices of the factorization scale /i = m, ^/Q' 2 + 4m 2 . The scale dependent terms in the NLL 



approximation to Eq. ( p0| ) have been multiplied with a damping factor 1/a/1 + rj in order to 
suppress spurious contributions from the large rj region, i.e. far above partonic threshold.^ 
Figs. |i~0| a andp^Ojb show that, for the chosen kinematical range in x and Q 2 , _P| harm as denned 



by Eq. (|50|) originates entirely from that region in r\ in which we approximate the exact results 
very well. In words, for the values of x shown, _P 2 = harm is completely determined by partonic 
processes close to threshold. We observe that the NLL approximation to Eq. ( |50"|) is much 
better than the LL one at the hadronic level as well. As we shall shortly show in more detail, 
it is also significantly more stable under variations of the factorization scale. 

In Figs. [TT| a and |Tl]b we display the results for Eq. (^) evaluated to NNLO in NLL approx- 
imation for the same kinematical values and choices for /i as before. Again, we have multiplied 
the NLO scale dependent terms c^'g and also all NNLO terms c^'g , / = 0, 1, 2 with a damping 



factor 1/^/1 + rj. Figs. 11a and lib reveal that the NNLO corrections are numerically quite 



sizable. However, the variations with respect to the mass factorization scale are still large for 
both values of x = 0.1, 0.01. 

Let us therefore turn to the issue of factorization scale dependence of the hadronic structure 



5 c 



See the discussion in section 3.2. If we add the damping factor also to the scale independent terms, our 



NLL results at fi = m are only 5% smaller. 
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Figure 11: (a): F% harm (x, Q, z max ) as a function of z max at NNLO in NLL approximation with the 
CTEQ4M gluon PDF, x = 0.1, m = 1.6 GeV, Q 2 = lOGeV and with the damping factor 1/y/T+Tj 
on all scale dependent terms and on c^' g . Plotted are: Scale choices \i = m (upper curve) and 
fj, = \/Q 2 + 4m 2 (lower curve), (b): Same as Fig. [TT| a for x = 0.01. Plotted are: Scale choices \x = m 
(lower curve) and \i = \J Q 2 + 4m 2 (upper curve). 



function in more detail. The scale dependence of _F 2 charm ; cf. Eq.(40), at NLO is exhibited in 
Figs. |l~2| a and |12|b over a range in m < fi < I^Q 1 + 4m 2 . For x = 0.1,0.01 we compare the 
exact result of Ref . |4]] with the LL and the NLL approximations from Eq. ([37]) . The coefficient 
c^'g is in the approximations again multiplied by the damping factor 1 / y/1 + r\. 

We see from Figs. 12a and |P2| b that resummation to NLL accuracy is crucial to obtain a 
reliable behaviour with respect to the /i-dependence of i^ harm at NLO. The NLL curve traces 
the exact result extremely well, in particular at phenomenological interesting smaller x values, 
when the support of gluon PDF extends to regions of tj — 1 to 10. In fact, Fig. [I2[ b clearly 
shows that in this kinematical domain, the LL approximation is rather poorly behaved. 

Continuing, we display in Figs. |l"3| a and [l3| b the effect of the NNLO corrections on the scale 
dependence of i^ 2 charm . We compare an improved NLL approximation to _P 2 charm at NNLO, which 
consists of the exact result of Ref. 0] for i^ 3 ™ at NLO and, additionally, our NLL approximate 
NNLO results from Eq . (|39|) , with the so-called best approximation to _F 2 charm at NNLO. This 
best approximation extends beyond the NLL accuracy at two-loop order as it contains the full 
exact results for c^'g and c[> 2 ^ from Eqs. (^6|) , (47), such that only remains approximate 
to NLL accuracy. 

We observe practically no improvement, with respect to the NLO result, in the scale depen- 
dence of _P| harm at large x near threshold (Fig. [T5|a), where the relative variations with respect 
to /j, were largest before. We also see here that the two approximations to _P 2 charm at NNLO, 
the NLL improved and the best, agree very well. They are differing only by terms beyond the 
NLL accuracy, which turn out to be rather small numerically. Therefore it is also unlikely that 
any corrections to c^f 1 beyond the NLL accuracy will help to soften the large //-dependence of 
F 2 charm at large x. 

On the other hand, for smaller x (Fig. [T~3|b) , we indeed see an indication of increased NNLO 
stability. The best NNLO approximation to _F 2 charm reduces the value of the relative variation 
F 2 charm (/i = m)/F 2 charm (/i = 2^Q 2 + 4m 2 ) from 1.13 at NLO to 1.07 at NNLO, so that there is 
practically no dependence on \x left. Naturally, for smaller x, where one probes regions away 
from threshold, terms beyond the NLL accuracy also have some numerical importance, as one 
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Figure 12: (a): The ^-dependence of the charm structure function i? 2 charm at NLO with the CTEQ4M 
gluon PDF, x = 0.1, m = 1.6 GeV and Q 2 = 10 GeV. Plotted are: The exact result (solid line), the 
LL approximation (dotted line) and the NLL approximation with the damping factor 1/^/1 + r/ on 
the scale dependent terms (dashed line), (b): Same as Fig. |l^ a for x = 0.01. 
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Figure 13: (a): The /U-dependence of the charm structure function i?| harm a t NNLO with the 
CTEQ4M gluon PDF, x = 0.1, m = 1.6 GeV, Q 2 = 10 GeV. Plotted are the improved NLL approxi- 
mation (dashed line; exact NLO result plus NLL approximate NNLO result with the damping factor 
1 /^/l + Vj) and the best approximation (solid line; exact results at NLO and for cJ> ' g 



NLL approximation c. 



and c 2 ' g plus 

2g with the damping factor + rj). (b): Same as Fig. 13a for x 
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,(2,2) 
"2,3 

0.01. 
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Figure 14: (a): The x-dependence of the charm structure function _p|harm a ^ jsjlo with the CTEQ4M 
gluon PDF, fi = m = 1.6 GeV and Q 2 = lOGeV. Plotted are: The exact result (solid line), the 
LL approximation (dotted line) and the NLL approximation (dashed line), (b): The x-dependence 
of the ratios F^ 0) /F^ } (solid line) and F^™ L0) /F^ 0) (dashed line) with if^ LO) in 
the improved NLL approximation (exact NLO result plus NLL approximate NNLO result with the 
damping factor 1/^/1 + rj) and parameter choices as in Fig. 14a. 




Figure 15: (a): The x-dependence of F| harm at NLO (lower band), scaled down by a factor of 
10 and (upper band) in the improved NLL approximation at NNLO (exact NLO result plus NLL 
approximate NNLO result with the damping factor 1/y/l + rj) with the CTEQ4M gluon PDF, fj, = m, 
Q 2 = 10 GeV. The bands correspond to the charm mass variation 1.35 GeV < m < 1.7 GeV. (b): 
Same as Fig. ICa for \i = \/ Q 2 + 4m 2 . 
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may see from the difference between the two approximations to i^ harm in Fig. |13|b. Finally, 
we note that a fully consistent NNLO scale dependence study would require not-yet-available 
NNLO parton densities, and the presently also unknown three-loop anomalous dimensions. 



The final comparison in this subsection involves the size and x dependence of _P| harm at a 
fixed value of the factorization scale fi — m. In Fig. [H| a we compare the exact result of Ref. 
with the LL and the NLL approximations from Eq. (|37|) , over a range of x, 0.003 < x < 0.3. 
We find that the deviations of the NLL approximation from the exact result are often very 
small, at most 10% for the value of x = 0.003. The excellence of the NLL approximation to 
the x dependence also holds for other values of \i. 

In Fig. |H| b we display for the same kinematics and over the same range in x the effect of 
the NNLO corrections. We plot the K-factors -F^awlo) /^2\nlo) an d> f° r comparison, also 
^2 C ^LO)/-^2 h (LO) 0- At NNLO we have taken the improved NLL approximation to F 2 charm ( we 
recall: the exact NLO result plus the NLL approximate NNLO result with the damping factor 
l/y/l + rj). We see that particularly for smaller x, the size of the NNLO corrections is negligible, 
the K-factor being close to one, whereas for larger x, their overall size is still quite big, almost 
a factor of 2 at x = 0.1. 

We note furthermore that the largest error on the approximate NNLO F 2 charm , as at NLO 



14|, is still due to the uncertainty in the charm mass. Varying the mass of the charm from 



1.35 to 1.7 GeV, we found sizeable variations in the absolute value of _F 2 charm ; in particular a 



larger x. In Figs. 15a and 15b we plot for two values of the factorization scale fx the improved 
NLL approximation to _P 2 charm as a function of x, 0.003 < x < 0.3. For comparison, Figs. 



and |i5[ b also contain the exact NLO result with same variations of the charm mass. In general, 
decreasing the charm mass increases i? 2 charin . 



charm 



This concludes our investigation of the inclusive structure function F 2 C 
3.4 dF£ haTia /dp T 

In this subsection we study the distribution <iF 2 diarm /<ipr- Soft gluon resummation might be 
especially fruitful here because the requirement that the detected charm quark has a fixed 
transverse momentum increases the sensitivity to threshold dynamics by effectively increasing 



the heavy quark mass from m to the transverse mass = ym 2 + p\. 

Eq. (|5|) provides us with the differential distribution d 2 F2/ ' dT\d\J\ from which we can derive 
the experimentally relevant distribution d 2 F2/dydpT, with y being the rapidity in the cm. 
frame of the vitual photon-proton system and pt the transverse momentum of the detected 
heavy quark. In terms of the Mandelstam variables of Eq. (Q), we can write the energy of the 
outgoing heavy-quark as 

E = _«!±^i. (B1) 

The transverse momentum pt and the longitundinal momentum pi are determined by 

S' 2 (p 2 T + m 2 ) = S'^ + Q^ + Q 2 ^, (52) 

Vl = E 2 -m 2 -p 2 T} (53) 

6 For ^^(LO) we use d a two-loop a s and NLO gluon density. 
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Figure 16: (a): The differential distribution dF^^/dpr as a function of px at NLO with the 
CTEQ4M gluon PDF, x = 0.01, m = 1.6 GeV, Q 2 



10 GeV and scale choice fi = y Q 2 + 4(m 2 + pj,). 
Plotted are: The exact result (solid line), the LL approximation (dotted line) and the NLL approxi- 
mation with the damping factor 1/^/1 + r\ on the scale dependent terms (dashed line), (b): Same as 



Fig. 16a for x = 0.001. 



where S' = S + Q 2 . The rapidity y 



2 \E-Pl, 



(54) 



may then be expressed in terms of (T 1; Ui) by Eqs. ( |51~D and (|53]) . The transformation (T 1; Ui) — > 
{pr,y) is then denned. Neglecting all light initial state quarks, one obtains from Eq. @ the 
result for the inclusive distribution dF^jdpx^ 



d 2 F 2 {x,p Tl Q 2 1 m 2 ^ 
dpi 



2 V+ 1 d 

I dy I ~^9/p( z ^ 2 ) ^2, ff (^,s 4 ,ti, Mi, Q 2 ,m 2 ,/i 2 ) , (55) 



where <p% < 5/4 — m 2 in the physical region and 



y 



± = ± cosh 



-if v 7 ^ 



2 JVTiT 



(56) 



With the above expressions, and using the one- and two-loop expansions of the hard part 
in section 3.1, we can construct NLL approximations to c?F 2 charm / 'dp^ at NLO and NNLO. In 
analogy to the inclusive ircharm^ we can ^ nen CO mpare the approximate NLO results with the 



exact ones of Ref. [[15], and, where appropiate, make NNLO estimates. 

In Figs. |16"| a and |16"l b we show the NLO results vs. px for x 
1.6 GeV and Q 2 



10 GeV. As in Ref. 15 



0.01 and 0.001, and m 
the factorization scale has been chosen as u 



Q 2 + 4(m 2 +p 2 ).Q 



At NLO, we compare our LL and NLL approximate results with the exact results of Ref. [|15 
and we see that the exact curves are reproduced well both in shape and magnitude by our NLL 



7 We have checked that the distribution dFZ h ' M ' 1 



^/dpx, when integrated over px, reproduces the inclusive 
charm structure function F^ haTnl i n section 3^ for the choice /i = \/ Q" 1 + 4m 2 . 
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Figure 17: (a): The differential distribution dF| tarm / 'dpT as a function of px with the CTEQ4M 
gluon PDF, x = 0.01, m = 1.6 GeV, Q 2 = lOGeV and scale choice y, = \JQ 2 + 4(m 2 + p T )- Plotted 
are: The exact NLO result (solid line) and the improved NLL approximation at NNLO (dashed line) 
(exact NLO result plus NLL approximate NNLO result with the damping factor 1/^/1 + rj). (b): Same 
as Fig. 0a for x = 0.001. 



approximations, whereas the curves for LL accuracy systematically underestimate the true 
result. In particular, we want to point out that for dF| harm /<ipT the kinematical region in 
which our NLL threshold approximation applies, extends well down to x = 0.001, with an error 
of at most 15% for small px, see Fig. |l6|b. The deviation at larger pt is presumably related to 
ln(pr/m) logarithms, associated with heavy quark fragmentation. 

For the same choice of kinematics, Figs. |l~7| a and [T7f b display the results for the improved 
NLL approximation to dF% harm /dp T at NNLO. We see that the small px region receives sizeable 
contributions, the value of the maximum increases by 40% - 50%, whereas in the large pt region 
there is little change. 



4 Conclusions 

We have performed the resummation of threshold logarithms for the electroproduction of heavy 
quarks, to next-to-leading logarithmic accuracy, and in single-particle inclusive and pair invari- 
ant mass kinematics. For the former, we have executed an extensive numerical investigation 
into the quality of the approximation for the inclusive charm structure function, and for its 
transverse momentum distribution. In addition we have provided both analytical and nu- 
merical results for NNLO approximations, for those kinematic configurations where the NLO 
approximation fits the exact results well. 

We found that the region of applicability of our analysis extends well into the kinematic 
range of the HERA experiments, a result of the reaction being largely driven by initial state 
gluons, which on average have not much energy. Our studies show a clear superiority of next-to- 
leading logarithmic threshold resummation over leading logarithmic resummation. In view of 
the present and much larger future HERA data sample, we hope that our results will find use, 
for example in the determination of the gluon density, and its uncertainty, and, more generally, 
that they provide motivation for further application and development of threshold resummation 
techniques. 
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Appendix A: Laplace transforms 

In this appendix we list the Laplace transform table needed to obtains the results of section 3. 
Define 

~\n n (s A /m 2 y 



o 



s 4 /m 2 

For the lowest four values of n this integral is, up to 0(1/N) 



(A.l) 



I (N) = -InN, (A.2) 
h(N) = Iln 2 A>+ic 2 , (A.3) 

h(N) = -iln 3 A>-C 2 miV- jjc 3 , (A.4) 
h(N) = iln 4 iV+^C2ln 2 iV + 2C3lniV+^C 2 2 + ^C4, (A.5) 
with N = Ne lE and 7^ denoting the Euler constant. Note that these answers are identical to 



the ones for the Mellin transform in 2£ . 



Appendix B: NLL resummed heavy quark pair-inclusive 

cross section 

Here we consider the singular behavior of the cross section for heavy quark electroproduction, 
in pair-invariant mass (PIM) kinematics. We present the resummed cross section, and its one- 
and two-loop expansions. Specifically, we consider the reaction 

nq)+P(p) — QW)+X( Px ), (B.l) 

We denote p' 2 = {jp\ + P2) 2 = M 2 and define r = M 2 /S. The definitions of other invariants are 
given in section 2. The cross section of interest is differential with respect to M 2 , the scattering 
angle 6 (in the pair center of mass frame), and the rapidity y of the pair. Adopting the same 
approximations as in section 2, it reads 

d 3 a(x,M 2 ,6,y,Q 2 ) 1 r r dx' , 2 / 1 1\ / xAm 2 /Q 2 \ 

dM 2 d cos 9 dy = ^J dZ J l^^ p{x ^ )6 { y -2 ln ^) 5 { Z -^rr^) 

x u(z,e^^,a s {p)) . (B.2) 



We have used here the observation |30| that the inclusive partonic cross section may be used 



to compute the singular behavior of the hadronic cross section at fixed (small) rapidity. 
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The kinematics of the reaction (|B.1|) near threshold is determined by the vector £ M = (1,0). 
(Recall that for 1PI kinematics = p^/m)- Near threshold it decomposes as 

» = (l-r) = ( 1 -')( T 1 i) + 2 ' ;S ' C 

1 



Wi 



1 — X 



w s - 



(B.3) 



Following the methods described in section 2, we find, in moment space 



ij(N,e,^,a,(M) 



H(e,^,a,(M)) S(l, a ,(f) 

M/N / 



x expj^^iV^^^M 2 )! expj j %2 Re {T s (e* s (//))} \ . (B.4) 

For simplicity, we have set the factorization scale \x = M, however it may be easily restored 
following the arguments leading to Eq. fl24|). The various functions in this expression are given 
explicitly in Eqs. (p5|)-(|27|) and ([3(]), in which 

Ul = -s'(l+/3cos0)/2, *i = -s'(l-/3cos0)/2, (B.5) 



with (3 = Jl — Am 2 / s. The lowest order hard part in Eq. ( |B.2| ), defined by 



^ (( V, M, 0) = s' 5(1 - r) c7 Bom (s', M, 0) , 
is related to its counterpart for 1PI kinematics in Eq. fl34j) by 



a BoTn (s', M, i 



2vr 2 «/? Born / 

-cr 9 „ | s , t\ 



s'Q 2 2 ' 9 V 2 v " 11 - 2 

The NLL approximation to the exact NLO correction with /i = M is 

cj (1) (s',0,M) ~ a Bor V, 0, M) , 
where 



(B.6) 

(l-/3cos0),«i = -^(l + /3cos0)^) . (B.7) 

(B.8) 



a a (M) 



7T 



2 a 



ln(l - 2) 
1 - ^ 

r2\ 



+ 



1 



1 - z 



C A -2 In 



1 — x 



+ ReL 



/3 



where 



V(l-z)~ 


= lim < 


1 - z 





-0(1-2-5) 



ln' +1 (5)5(l-^) 



Z + l 



The NNLO corrections with ijl — M are to NLL accuracy 

J 2 \s',6,M) ~ Lf( 2 > a Born (s', 0, M) 

where 

a 2 (M) 



(B.9) 



(B.10) 



(B.11) 



7T 



-6 In 



2C 2 



1 — x 



"ln 3 (l - z) 




"ln 2 (l-^)" 




+ 


1 - z 




1 - z 


+ 



+ ReL 



P 



3 In 



'M 2 



rrr 



-2C A (& 2 + 3C F (ReL /3 + l)) 



(B.12) 
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